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Abstract: In this paper, we study the global existence of steady subsonic Euler flows 
through infinitely long nozzles without the assumption of irrotationality. It is shown that 
when the variation of Bernoulli's function in the upstream is sufficiently small and mass 
flux is in a suitable regime with an upper critical value, then there exists a unique global 
subsonic solution in a suitable class for a general variable nozzle. One of the main difficulties 
for the general steady Euler flows, the governing equations are a mixed elliptic- hyperbolic 
system even for uniformly subsonic flows. A key point in our theory is to use a stream 
function formulation for compressible Euler equations. By this formulation, Euler equations 
are equivalent to a quasilinear second order equation for a stream function so that the 
hyperbolicity of the particle path is already involved. The existence of solution to the 
boundary value problem for stream function is obtained with the help of the estimate 
for elliptic equation of two variables. The asymptotic behavior for the stream function is 
obtained via a blow up argument and energy estimates. This asymptotic behavior, together 
with some refined estimates on the stream function, yields the consistency of the stream 
function formulation and thus the original Euler equations. 
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1 Introduction and Main Results 



Multidimensional gas flows give rise many outstanding challenging problems. Since the 
solutions for the unsteady compressible Euler equations develop singularities in general [29]. 
it is not yet known which function space is suitable to study their wellposedness[2^. It 
is natural to start from the steady Euler equations to understand some important true 
multidimensional flow patterns. However, the steady Euler equations themselves are not 
easy to tackle, since the equations may not only be hyperbolic or hyperbolic-elliptic cou- 
pled system, but also have discontinuous solutions such as shock waves and vortex sheets. 
Therefore, a lot of approximate models were proposed to study fluid flows. An important 
approximate model is the potential flow, which originates from the study for flows without 
vorticity. Since 1950's, tremendous progress has been made on the study for potential flows. 
Subsonic potential flows around a body were studied extensively by Shiffman[28j, Bers[2l[3], 
Finn, Gilbargp^ I14j. and Dong[TO], et al. Subsonic-sonic flows around a body were estab- 
lished recently by Chen, et al [5] via compensated compactness method. Signiflcant progress 
on transonic flows was made by Morawetz. She flrst showed the nonexistence of smooth 
transonic flows in general [2T | 1221 123^ I24j . and later worked on existence of weak solutions to 
transonic flows by the theory of compensated compactness |25|, I26j. Existence and stability 
of transonic shocks in a nozzle for potential flows were achieved recently with prescribed 
potential at downstream in [6l l^. Xin and Yin obtained existence and nonexistence of 
transonic shocks in a bounded nozzle with prescribed pressure at downstream was obtained 
in |34t l35j. Recently, well-posedness for subsonic and subsonic-sonic potential flows through 
inflnitely long 2-D and 3-D axially symmetric nozzles, was established in |3H I32j. For the 
study on other aspects on subsonic potential flows, please refer to [T2 l [T5 l [T6]. 

Besides the potential flow, there is another important approximate model to compress- 
ible Euler equations, incompressible Euler equations, which approximate to the compress- 
ible Euler equations for flows with small Mach numbers. For the study on the existence of 
steady incompressible Euler flows in a bounded domain, please refer to (UEOjIIB], etc, and 
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references therein. 

For the full compressible Euler equations, well-posedness and nonexistence of a transonic 
shock in bounded nozzles with prescribed pressure at downstream has been obtained in 
[33\ [36l [T9l [20]. Existence and stability of transonic shocks in nozzles with prescribed 
velocity at downstream was shown in [3]. 

In this paper, we study the existence of global steady subsonic Euler flows through 
general infinitely long nozzles. 

Consider the 2-D steady isentropic Euler equations 

{pu)x, + {pv)x2 = 0, (1) 
{pu^)xi + ipuv)x2 +Pxi = 0, (2) 
{puv)x^ + {pv^)^.-^ + Px^ = 0, (3) 

where p, {u,v), and p = p{p) denote the density, velocity and pressure respectively. In 
general, it is assumed that p'{p) > for p > and p"{p) > 0, where c{p) = \/p'{p) is called 
the sound speed. The most important examples include polytropic gases and isothermal 
gases. For polytropic gases, p = Ap^ where ^ is a constant and 7 is the adiabatic constant 
with 7 > 1; and for isothermal gases, p = c^p with constant sound speed c [U]. 
We consider flows through an infinitely long nozzle given by 

Q. = {(2;i,X2)|/i(2;i) <X2 < f2{xi), -00 < xi < 00}, 

which is bounded by Si = {{xi, X2)\x2 = fi{xi), —00 < xi < 00}, (i = 1,2). Suppose that 
Si{i = 1,2) satisfy 

/2(a:i) > /i(xi) for xi £ (-00, cx)), (4) 
/i(xi)^0, /2(a:i) ^ 1, as xi ^ -03, (5) 
/i(xi)^a, /2(xi)^6>a, as xi ^ +00, (6) 

and 

||/i||(72,Q(]K) < C for some a > and C > 0. (7) 
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It follows that satisfies the uniform exterior sphere condition with some uniform radius 
r > 0. 

Suppose that the nozzle walls are impermeable solid walls so that the flow satisfies the 
no flow boundary condition 

{u,v) ■ n = on 50, (8) 
where n is the unit outward normal to the nozzle wall. It follows from ([T|) and ([8]) that 

{pu, pv) ■ ndl = m (9) 



I 

holds for some constant m, which is called the mass flux, where I is any curve transversal 
to the direction, and n is the normal of / in the positive xi-axis direction. 

Due to the continuity equation, when the flow is away from the vacuum, the momentum 
equations are equivalent to 

UUx^ + VUx2 + Kp)xi = 0, (10) 
nuxi + vvx2 + h{p)x2 = 0, (11) 

where h{p) is the enthalpy of the flow satisfying h'{p) = p'{p)/p- So h{p) is determined up 
to a constant. In this paper, for example, we always choose h{0) = for polytropic gases 
and h{l) = for isothermal gases. After determining this integral constant, we denote 
Bo = infp>o h{p). 

It follows from ([10]) and ([n]) that 

{u,v)-V{h{p) + ^{u^ + v')) = 0. (12) 

This implies that "^^^^ + h(p), which will be called Bernoulli's function, is a constant along 
each streamline. For Euler flows in the nozzle, we assume that in the upstream, Bernoulli's 
function is given, i.e., 

1/ 2 -I- 71^ 

+ h{p) — > B{x2) as xi — oo, (13) 
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where B{x2) is a function defined on [0, 1]. 

Now let us state our main results in the paper 



Theorem 1 Let the nozzle satisfy and B_> Bq. There exists a 5q> {) such that if 

inf B{x2) = B, B'{0) < 0, B'{1) > and ||S'(x2)||co,i([o il) = < ^^o, (14) 

a;2e[0,i] 

then there exists m > 25^^ such that for any m £ ((5^/^,m), 

1. (Existence) there exists a flow satisfying the Euler equations ([2P-(0j; the boundary 
condition mass flux condition (0), and the asymptotic condition M3^) : 

2. (Subsonic flows and positivity of horizontal velocity) the flow is globally uniformly 
subsonic and has positive horizontal velocity in the whole nozzle, i.e., 

sup(?x^ + f ^ — (?{p)) < and u > in Cl; (15) 
n 

3. (Regularity and far fields behavior) Furthermore, the flow satisfies 

llp|lci'°^(n)> ll^llci-«(n)) lbllci.°(n) < C (16) 
for some constant C > 0, and the following asymptotic behavior in far fields 

/? — > po > 0, u ^ uo{x2) > 0, w — > OS xi — > — oo, (17) 
Vp ^0, V-u ^ (0,Uo(x2)), Vv ^ as xi ^ -oo, (18) 

uniformly for X2 G Ki (s (0, 1), and 

p Pi > 0, u ^ Ui{x2) > 0, V — > OS xi ^ +00, (19) 
Vp ^0, V-u ^ (0,u'i(x2)), V^; ^ as xi ^ +oo, (20) 

uniformly for X2 € K2 (e (a, 6), where po and pi are both positive constants, and po, 
pi, uq, and ui can be determined by m, B{x2) and b — a uniquely; 

4- (Uniqueness) the Euler flow which satisfies (Cp-f^), boundary condition (0), asymp- 
totic condition ^3\} . mass flux condition (0), U5\) . and asymptotic behavior 1 l'/)-[T3\) 
is unique; 
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5. ( Critical mass flux) If, besides B also satisfies 

B'{0) = B'{1) = 0, (21) 

then rh is the upper critical mass flux for the existence of subsonic flow in the following 
sense: either 

sup(u^ + v'^ — (?{p)) as m ^ rh, (22) 

n 

or there is no a > such that for all m S (rh, rh + a), there are Euler flows satisfying 
(l2P-(0j; boundary condition asymptotic condition mass flux condition 

^^). and asymptotic behavior (|T7|j-( l7^] and 

sup sup(c2(/9) - (n^ + v"^)) > 0. (23) 

There are a few remarks in order: 

Remark 1 Here we obtained only the existence of the Euler flows in the nozzle, and the 
uniqueness in a special class of flows, but not the uniqueness for general Euler flows. For 
the issue on the uniqueness for steady incompressible Euler flows in a bounded domain, 
please refer to f30| . 

Remark 2 It can be shown by modifying the analysis in this paper slightly without further 
difficulties that there exists a subsonic full compressible Euler flow in the nozzle, if the 
entropy is prescribed in the upstream. 

Remark 3 The subsonic Euler flows in half plane was studied in [8j recently. Although 
stream function formulation is also introduced in [8j, however, the far flelds conditions are 
different from ours. Furthermore, we obtain critical upper bound of mass flux for existence 
of subsonic flows in nozzles. 

The rest of the paper is arranged as follows: in Section [3l we reformulate the problem 
by deriving the governing equation and boundary conditions for Euler flows in terms of a 
stream function, provided that the Euler flow has simple topological structure and satisfies 
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the asymptotic behavior (fT7]) - p^ . In Section [3l existence of solutions to a modified elhptic 
problem is established. Subsequently, in Section [H we will study asymptotic behavior 
of solutions in a larger class and show uniqueness of the solution to the boundary value 
problem. The existence of boundary value problem for the stream functions will be a 
direct consequence of these asymptotic behavior and uniqueness. In Section [5l some refined 
estimates for the stream function will be derived. Combining these estimates with the 
asymptotic behavior obtained in Section [4] will yield the existence of Euler flows which 
satisfy all properties in Theorem [TJ Finally, in Section [H we will show the existence of the 
critical mass flux. 



2 Stream-Function Formulation of the Problem 

We start with some basic structures of the steady Euler system. The steady Euler system 
([T])-(l3]) can be written in the following form, 

AU,, + BU,, = 0, 



where 



A 



( uc^p) 

p 



c\p) 



\ 



V 

Let A be the solution of 



pu 
puj 



B 



p 




c^{p)^ 



pv 



\c^{p) pv J 
det{XA -B) = 0. 



U 



It follows from straightforward computations that (j24p has three eigenvalues 



(24) 



V uv ^ c{p) ^/v?+^P^^^(P\p) 

Ai = — , A± = 5 5 . 

Therefore, at the points where + v"^ — (?{p) > 0, i.e., the flow is supersonic, (j24p has 3 
real eigenvalues, the Euler system is hyperbolic. When u'^ + v"^ — (P'{p) < 0, i,e., the flow 
is subsonic, (j24p has a real eigenvalue and two complex eigenvalues, the Euler system is a 
hyperbolic-elliptic coupled system. Therefore, even for globally subsonic flows, one has to 



resolve a hyperbolic mode. Moreover, for flows in infinitely long nozzles with both ends at 
infinity, it seems difficult to get uniform estimates for hyperbolic mode. 

To overcome the difficulties mentioned above, we introduce the stream functions for the 
2-D steady compressible Euler flows, and derive an equivalent formulation for Euler flows 
in terms of the stream functions when the flow satisfies certain asymptotic behavior. 

It follows from ([10]) and (Ull) that 



9x2(^x1 + vux2 + Hp)xi) - dx^{uvx^ + vvx2 + h{p)x2) = 0, (25) 
which yields that 

(n, v)-\luj + a;div(n, v) = 0, (26) 
where oj = Vx^ — is the vorticity of the flow. By the continuity equation ([T]), one has 

(u, v) ■ Vuj + wdiv(ti, v) = (u, v) ■ (Vlv — ^^^ ) = p(u, v) ■ V ( — ) . 

P VP/ 

Therefore, away from vacuum, (j26p is equivalent to 

(n,z;)-vQ = 0. (27) 

We have the following proposition. 

Proposition 2 For a smooth flow away from vacuum in the nozzle satisfying ^ and 
(0), the system consisting of ([T]), and ^27^ is equivalent to the original Euler equations 
([2P-(0i, if the given flow satisfies no flow boundary condition, 

u>0 inn, (28) 

and the following asymptotic behavior 

u, p and Vx2 <ire bounded, while v, Vx-^ and px2 0, as xi —oo. (29) 

Proof: From previous analysis, it is easy to see that smooth solutions to the Euler equations 
(d])-© satisfy {I]), {E]) and (l27|). On the other hand, it follows from (P), dS?]) and the above 
derivation that ()25p holds. Therefore, there exists a function <I> such that 



^xi = uuxi + vux2 + h{p)xi , = uvxi + vvx2 + h{p) 



X2- 



So, ([T2]) is equivalent to 

{u,v)-V<^ = 0. (30) 

Due to the no flow boundary condition ([81), $ is a constant along each component of 
the nozzle boundary. If, in addition, 

— > as xi — > — oo, (31) 

then $ — > C as j;i — > — oo. On the other hand, it follows from (j28|) that through each point 
in O, there is one and only one streamline satisfying 

< 

V 

which can be defined globally in the nozzle (i.e., from the entry to the exit). Furthermore, 
it follows from ([T]) that any streamline through some point in 0, can not touch the nozzle 
wall. Suppose not, let the streamline through (x5,X2) pass through (xi,/i(xi)). Due to ([I]) 
and no flow boundary condition, one has 

/ {pu){Xi,X2)dX2 = 0. 

This contradicts ([28]) . 

Thus, one can always solve (j30|) in the whole domain 0, which yields 

$ = C in 17, 

if dSID holds. Therefore, = = in 0, i.e., ([lO]) and ([11]) hold globally in the nozzle. 
Thus, both ([2]) and ([3]) are true. It is obvious that holds if ([29]) is vahd. □ 

It suffices to prove the existence of solutions to the system ([T|), (fT2]) and ([271) satisfying 
(|28[) and (|29D . 

However, system ([TJ, (I12p and (I27p is not easy to study directly either, since for in- 
finitely long nozzles with both ends at infinity, it seems difficult to estimate the solutions 
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for transport equations (jl2p and (j27p . Instead, we will use an equivalent formulation for 
©, ([H and (IIZD. 

It follows from ([1]) that there exists a stream function ^ such that 

Thus for the flows away from the vacuum, (f27|l is equivalent to 

V^V' • ^ (^) = 0' (32) 

where V-^ = (-S^^a.^^^J- Note that ([32]) means that — and ij: are functionally dependent, 
therefore, one may regard ^ as a function of ip. Set 

- = WN)). (33) 
P 

Similarly, p2p is equivalent to 

therefore, h[p) + \ {u^ + w^) is also a function of ^/^. We define this function by 

h{p) + \{u^ + v^)=B{i:). (34) 

Furthermore, it follows from ([8]) that the nozzle walls are streamlines, so "0 is constant on 
each nozzle wall. Due to ([9]), one can assume that 

V' = on 5*1, and = m on 82- (35) 

In order to get an explicit form of B, first we study the density-speed relation via 
Bernoulli's law ()34|) carefully. 

Note that p'{p) > for p > and p"{p) > 0, therefore h'{p) = p'{p)/p > for p > and 
for some fixed p > 0, 

h{p)=h{p)+ r?^ds>h{p)+ r ^^dsioT p> p. 

Jp S Jp s 
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This yields that h{p) — oo as p — oo. On the other hand, since infp>o h{p) = Bq, h{p) — Bq 
as p — 0. Thus for any s > Bq, there exists a unique g = g{s) > such that 

h{e{s)) = s. 

Moreover, for the state with given BernouUi's constant s, the density and speed satisfy the 
relation. 

Hp) + \ = s. 

Therefore, the speed q satisfies 

q = ,/2{s-hip)). 

Hence, for fixed s, g is a strictly decreasing function of p on [0, g{s)]. By the definition of 
g{s), one has q{g{s)) = < c{g{s)). Now we claim that q{0) > c(0). Indeed, one can prove 
this claim in two cases. First, if c(0) > 0, then 

hip) = hip) + / ^—-^ds < hip) + / -^ds ^ -oo as p ^ 0. 

Jp S Jp s 

Therefore, q{0) oo. Thus q{0) > c(0). Second, if c(0) = 0, q{p) ^/2(F^^Bo) as p ^ 0, 
therefore, q{0) > = c(0). This completes the proof of the claim. Since c^(p) =p'{p) is an 
increasing function of p, there exists a unique g{s) G [0, g{s)] such that 

c'{g{s)) = q'{g{s)). 

In summary, for any given s > Bq, there exist g = g{s), g = g{s) and T = T{s) such that 

h{g{s)) = s, h{g{s)) + = s, and c\g{s)) = F^s), (36) 

where g{s), g{s), and T{s) are the maximum density, the critical density, and the critical 
speed, respectively for the states with given Bernoulli's constant s. Set 

E(s) = g{s)^/2is-h{g{s))). (37) 

Then direct calculations show that 

dg g dg 1 
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and 



Thus 



dT. _ y^2{s - h{g{s))) 1 2p' {Q)+fp" (g) 



> 0, ^ > 0, and — > 0. 
as ds ds 

Obviously, g{s) < q{s), if s > Bq. By the continuity and monotonicity of q(s) and g(s) 
there exists a unique 5>0 such that 



giB + 5) = giB)- 



(38) 



(39) 



Moreover, it follows from (j36p that there exists a uniform constant C > such that 

C'^ < q{B) < g{B) = g{B + S) < g{B + S)<C, 
C-^ < q'{s) < C, C-^ < g'{s) <C, if s G {B,B + 5), 

< h'ip) <C,iipe ig{B), m + 6)), 

< S(s) < C7, if s e {B,B + 5). 

Later on, C will to denote generic constants which depend only on B and ^, and thus 
essentially on B^. 

In order to study the relationship between density and mass flux with given Bernoulli's 
constant, let us investigate the function defined by 

I{p) = 2p\s-h{p)). 

Direct calculations show 



dp 



4p(, - h{p) -p'{p)/2) = 2p{q\p) - c\p)). 



Therefore, for p G (0,£»(s)), ^ > 0; and ^ < for p G {q{s) , g{s)) . Moreover, /(O) = 
I{q{s)) = 0- Thus I{p) > if /9 G (0, g{s)) and / achieves its maximum at p = g{s). So, for 
fixed s, the relation 



Kp) + 



2p2 



(40) 
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defines a function Ai = I{p) wliich attains its maximum M = at p = g{s). Thus p IS 

a two- valued function of for 7W G [0, 1]^(s)). Denote the subsonic branch by 

p = J{M) for M £ (0,s2(s)), 

which satisfies J{M) > q{s). When s varies, this branch will be denoted by 

p = J{M, s) for {M,s)e {{M, s)\M G (0, ^^(s)), s > Bq}. (41) 

To determine the explicit form of W and B, one may study W and B in the far fields of 
the nozzle where the flow may have certain simple asymptotic structure. Indeed, for flows 
satisfying the asymptotic behavior (fT7|) - ([2U|) . one can determine po, pi, uo(x2) and ui{x2) 
first. Suppose that the flow satisfies (fT7|) . Then 

Kpo) + = B{X2), uo{x2) > 0, (42) 



and 



hold, which shows that 



and 



If B(x2) satisfies 



then 



1 

Pouo{x2)dx2 = m (43) 







uo{x2) = ^/2{B{X2) - h{po)), (44) 



m= / poy/2iB{x2) - h{po))dx2. (45) 







mf B{X2)=B, ||i?'(x2)|bo.i([o,i]) < 6, (46) 

3:2G[0,1] 



B= sup B{x2) <B + 5. (47) 

X26[0,l] 

Let S < 6/2. Then it follows from ([38]) that g{B{x2)) < q{B) < g{B). To obtain a global 
subsonic flow in the nozzle, it is necessary to show that for given B2{x2) and ra, (j45p has a 
solution satisfying pQ G (giB) , q{B_)) . Direct calculations yield that 

poy/2iB{x2) - hipo))dx2 < 0, for po G {giB),Q{B)). 
13 



dpo Jo 



It follows from ([39]) and dHj) that 

giB)y/2iB{x2) - hiQiB)))dx2 = C q{B)V'^{B{x2) - B)dx2 < C5^/\ 
Jo Jo 



In addition, 

'1 



g{B)^2iBix2)-hiQ{B)))dx2 
> j\{B)^j2{B-h{Q{B)))dx2 



Q{B)^2{h{Q{B)) - h{Q{B)))dx2 

1 



Q{B)p{h{Q{B + S)) - h{Q{B)))dX2 
> [ q{B) y/2 {h{Q{B + 5))- h{Q{B + 6/2)))dX2 

Jo 

Therefore, for any 7 G (0,1/3), there exists (5o G (0,5/2) such that (|15]) admits a unique 
solution po £ (oiB), ^(^S)), provided that < S < 5o and m G {6"', mi), where mi satisfies 
mi > C^^6}^'^ > 26q^'^ . Later on, for definiteness, we will choose 7 = 1/4. However, all 
results hold for 7 G (0, 1/3). 
By virtue of (|45|) . one has 



m 



po ^/2{B{x2) - h{po))dx2 

= f Po^/2{B{x2) -B + B- h{po))dx2 
Jo 

< C I ^/2{5 + h{Q{B))-h{pQ))dx2. 
Jo 



Thus 



5 + h{Q{B)) - h{po) > C-^m^ > 
Note that 7 < 1/3, therefore, there exists 60 G (0, 5o) such that if < 6 < do, then 

h{Q{B)) - Hpo) > C-H^\ 
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Consequently, if ||-B'(a;2)||(70,i(jo,i]) = S < 60, by virtue of ([39|) . there is a positive constant 
C such that 

' < g^B) -P0<C, 

C~^S^ <uo<C, 



1^^0(3^2)1 



B'{X2) 



(48) 



^2{B{x2)-h{po)) 

^ Wo{^2)]co^m])<c{6^~'^ + s^--'^). 

To determine the states in the downstream, we parametrize the streamhnes in the down- 
stream by their positions in the upstream. Due to p!7|) . (fT9]) . and (p8|) . we can define 



such that 



y = yis) for s G [0, 1] 



Hpo) + = Mpi) + ^i(y(^)) > 0, 



PQUo{t)dt 



piui(t)dt, 



(49) 

(50) 

(51) 
(52) 



y{0) = a, y(l) = b. 

The meaning of y{s) is that the streamUne which starts from (—00, s) will flow to (00, y{s)) 
The map (|49|) is well-defined since (|28p ensures a simple topological structure of streamlines 
It follows from (j5ip that 

Po^^o(s) = Piui{y{s))y'{s). 



Hence, 



ds 



Po"o(s) 



piV2(h(po)+#^-MPi)) 



(53) 



(54) 



where the parameter pi satisfies 



y(0) = a, 

Pouois) 



zds = b — a. 



PiV2(%o) + ^-%i)) 
It remains to show that there exists a pi £ {q{B), g{B_)) satisfying (|55p . By direct calcula- 
tions, one has 



(55) 



_d_ 
dpi Jo 



Pouois) 



PiV2(%o) + ^-%i)) 



ds > 0, for pi e {Q{B),g{B)). 
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First, there exists 5o £ (0,(5o) such that 

Pouois) 



f 

Jo 



-q{B - 6)\/2ih{po) + ^ - hi-g{B - 6))) 
Pouo{s) 



zds 







eiB - 5)y/2{B{s) -B + hiom - h{0{B - 5))) 



zds 



provided 5 < 5q. On the other hand, 



zds 



e{B)^2{h{p,) + ^-^-h{Q{B))) 
m 



< 



Q{B)^2{B-h{e{B))) 
m 

Q{B)^2{h{g{B)) - h{e{B))) 
m 



Q{B)J2[h{e{B + 5))-h{Q{B))) 



m 

< 



< 



q{B)^2 {h{e{B + 6)) - h{Q{B + 5/2))) 
m 



So there exists a unique p\ G {q{B_) , g{B_)) such that ([55]) holds, if < 6 < 6o and m € 
((^'^,m2) for some m2 > mm{mi, C'^^ {b — a)6}^'^} . Furthermore, one can choose 6o smaher if 
necessary such that m2 > 25j/^ As soon as pi is determined, y{s) and ui can be obtained 



from (IMl) and (fSOl) . 

Let us summarize the above calculations in the following proposition: 

Proposition 3 Let S > Bq. There exists 6o > such that for any B G C-'^'^([0, 1]) 

_ ~l/8 

satisfying with 6 < 6q, there exists fh > 26q such that 

1. there exist solutions (po, uq) to [4^-(43^ o.i^d {pi,ui) solving ^5U^-^5^ ifmG 

(51/4 

2. PO, pi G {q{B),q{B)); 

3. {po,uo) satisfies (4^! 

4- either po q{B) or pi q{B) as m ^ fh; 
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where B = max^.2g[o^i] B{x2). 

Proof: We need only to verify the last statement. 

First, if m G {5^^^ ,m2)^ both po ^-nd pi belong to {q{B) , q{B_)) . For a given B{x2), as 
m increases, po decreases. If m — > m = Jq^ q{B)\J2{B{x2) — h{g{B)))dx2, then po — ^ q{B). 
Therefore there exists an upper bound for m to ensure the existence of po, pi G (oiB), q{B_))- 
Define 

fh = sup{s|m S {5'^,s), there exist po,Pi G (^/(-S), ^(5))} (56) 

Obviously, fh £ [1712, m]. Note that po and pi are uniformly away from g{B_). If neither po 
nor pi approaches to g{B) as m ^ fh, then there always exist po, pi G {q{B),q{B_)) for 
m G {5^,fh + e) with some small e > 0. This contradicts with the definition of m. So the 
proof of the Proposition is finished. □ 

We now can determine W and B in the upstream. Suppose that the flow satisfies the 
asymptotic behavior (I17p . Then in the upstream, a stream function can be chosen so that 



X2 

^ = I PoUo{s)ds, 
10 



(57) 



and < < m. Since pouo(s) > for s G [0, 1], '(/' is an increasing function of X2. Thus 
one can represent X2 as a function of ^, 



Define 



X2 = k(V'), < V < m. 



fiij) = n'o(K(V)), and F(V') = txo(K(V)). 



(58) 



Then / and F are well-defined on [0,m]. 




The parametrization of flows by a stream function 
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It follows from the proof of Proposition [2] that through each point (xi,X2) G there 
is one and only one streamline which starts from the entry, provided that ()28p holds in fi. 
By the definition of streamlines, along each streamline, the stream function is a constant, 
therefore, through any (xi,X2) in the nozzle, there exists a unique streamline originated 
from {—oo, with ■0 = ■0(xi,j;2)- Since Bernoulli's function is also invariant along a 

streamline. 



(Mp) + ^)(^i,^2) = (Mp) + 



Thus in the nozzle, one has 



)(-oo,k(V')) = h{po) + 



F^^{xi,X2)) 



Similarly, by virtue of (j33[) . one has 



n{p, |V0|2, 0) = h{p) + ^ - /i(/,o) - = 0. 



/(V^(a;i,a:2)) 
Po 



P 



provided that (j28p holds. Furthermore, note that (j28|) implies 



(59) 



(60) 



< ip < m. 



(61) 



Thus both (|59|) and (j60j) do make sense. 

Next, we study the relationship between F and /. In the upstream, 

Tp = pouo{s)ds, 
Jo 



which yields 



So, ([58|) shows 



this implies 



k'(V) 



1 



pouoiKii^)) PoF{iIj) 



F'{i;)=u'^{KmK'{^)=f{i^) 



1 



poF{i,) 



f{i;)=poFi^)F'i^). 



(62) 
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Furthermore, if B satisfies (|14p with < 6 < Sq, and m £ (5'^ ,m), then 



0-^6"/ <F<C, 

F'{m) > and F'{0) < 0, 



"o(«(^)) 



\F'm 

[^'(V')]c0.l([o,m]) 

It follows from ()4ip and (|59p that the subsonic flows in the nozzle satisfy 

p=J{\V^P\^h{po) + ^^), 



(63) 



(64) 



if they have asymptotic behavior ()17p . Furthermore, by the definitions of vorticity and 
stream function, one has 



CO 



-div 



P 



Thus, the stream function satisfies 



div 



where 



if(|VV^P,^) = J(|VVp,%o) + 



(65) 



(66) 



Our major task in the rest of the paper is to show the existence of solutions to the 
following boundary value problem 



div 



j^^^l^j =F(V^)F'(V^)i/(|VV^|2,V^) in 



X2-fl{xi) 



(67) 



/2(a:i)-/i(a;i) 

and show that the flow field induced by 



m on 517, 



P P 



satisfies ()17p -(|20 p . We will obtain further the estimates ()6ip and ()28p for the solution to 
(f67ll in order to get the existence of Euler flows. 
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Remark 4 It is easy to see that if B = B_, i.e., the flow has uniform Bernoulli's constant, 
then, F is a constant and the equation ([65|) reduces to 



This is nothing but the potential equation. Therefore, it is reasonable to use (j67p to formu- 
late the problem for the Euler flows through the nozzles as a perturbation of the potential 
flows. The condition m > (5'^ (7 < 1/3) ensures 



which guarantees that the magnitude of vorticity \FF'\ is sufficiently small, thus one can 
regard the potential flow as a leading ansatz for the Euler flow. 

3 Existence of a Modified Boundary Value Problem for Stream 
Function 

There are two main difficulties to solve the problem (j67p . The first difficulty is that the 
equation in (j67p may become degenerate elliptic at sonic states. In addition, H is not 
well-defined for arbitrary ip and IVV'I- The second difficulty is that this is a problem in an 
unbounded domain. Our basic strategy is that we extend the definition of F appropriately, 
truncate {Vipl appeared in in a suitable way, and use a sequence of problems on bounded 
domains to approximate the orginal problem. In this section we first get the existence of 
a modified problem on the unbounded domain, which indeed solves the original problem 
together with the asymptotic behavior established in the next section. 




F\ > C-^S^ and < 6^'^^, 



Set 



if < s < m 



F'{m)(2m — s)/m, if m < s < 2m, 



= < 



F' 



(0)(s + m)/m, 



if — m < s < 



0, 



if 5 > 2m, or s < — m. 



V 
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It is obvious that g G C°'^(M) and 



ll5(s)llcO(Ri) < ll^'('S)||cn([OH)- 



Moreover, it follows from (j63p that 



g{s) > if s > m and g{s) < if s < 0. 

Furthermore, it follows from ||i^'(s)||cO([o,m]) ^ Cd'^^'^'^', ||-F"('S)||cO'i([o,m]) ^ CS'^^^'^' and 
m > 6^, that 

rg{s)\\co.i^^i)<C5'-'^. (69) 

Define 



F{s) = F{0) + / ~g{t)dt. 
Jo 

Then F' = g and F G C^'^(M). Moreover, because m > 6'^, there exists a suitably small Si 
such that when 6 < 6i, 



Bo<B-eo< h{po) + <B + eo 



(70) 



holds for some Eq > 0, where = sup^^ep,!] B{x2)- Moreover, (p3]l and ([Mil imply 



|F'lbo(Mi) < C5'-^^ and ||F'|bo.i(Mi) < CJ^"'^. (71) 



In the rest of the paper, we will always use the following notations 

^i(|vvP,V') 



^^(lvv^l^v), H.imM.) = ^(ivvp,^). 



It follows from direct calculations that 



i7i(|VV'|',V') 



.2 IVt/'P^ 



2p(c- - 



may go to negative infinity when the flow approaches sonic from subsonic. 
Choose a smooth increasing function (^o such that 



Cois) 



s, if s < — 2eo, 
-eo, if s > -£o- 
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Then define 

A(|VV^|2, V') = CodVVp - ^\B{m + S'(^W), (72) 
where S is the function defined in (1371) and 



Set 



B{i^) = h{p,) + —^. (73) 



mvi^Wi^) = J(A(|vvl^ v),%o) + ^^), (74) 



2 

where J is the function defined in (1411). A direct calculation shows 



^i(|vvP,V)- 



2(ii'2c2 _ A) 

Obviously, there exist two positive constants A(eo) and A(eo) such that 

<A,,{q,z)i,i,<m^ (75) 
holds for any z G R-^, g G and ^ G M^, where 

4-(g, z) = ^(|g|2, z)5i, - 2#i(|g|2, z)qiqj. (76) 
Instead of ()67p . we first solve the following problem 



(77) 



div(^^^l^) =F(V^)F'(V^)F(|VV'P,V') in 

^ /2(2:i)-/l(a::i) 

Proposition 4 iei i/ie boundary offl satisfy (0i-(0j- T/ien i/iere exists < Si < mm{5o,5i} , 



where 5o is defined in Section\^ such that if ||i?'||(70,i([o,i]) = (5 < 5i and m £ {6"', mi) with 
mi = 28^^"^ < m, where m is defined in / f56]) in Section then the problem |771 ) has a 
solution ip G C'^''^{^) satisfying 

< C(eo, S), IVV-P < T.'^iR - eo) - 2eo for some Eq > 0. (78) 

Proof: Note that the equation (j77p is uniformly elliptic and the domain is unbounded, 
one can use a sequence of boundary value problems on bounded domains to approximate 
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it. The key point is to obtain the estimate (j78p . Therefore, we first solve the fohowing 
boundary value problem 

^-(^c^) =^(^)^'(^)^(l^^l''^) (79) 

where VLl satisfies {(a^i, j;2)|(2;i, X2) G fi, — L < xi < L} C I^l C {(a^i, 2:2)1(2^15 2:^2) S 
— 4L < xi < 4L} for VL G N. Furthermore, one may choose Q.^ so that fi^, G C^'"^(0 < 

ai < a) satisfies the uniform exterior sphere condition with uniform radius rg, < ro < r, 

for all L > Lq with some Lq sufficiently large. For the explicit construction of such i^L, 

please refer to Appendix in [31| . 

The equation in ()77p can be rewritten as 

4-(^^,V')9.,V-^2(|VVp,V)|VV'|2 = F(V')F'(V')F=^(|V^|2,V'), (80) 



where 



^2(1 WP,^) = l(^)^'(^)H{\Vij\\,p){H' + SS'(C^ - 1)) 



i72(|v^|2,V;)c2- A(|VV'p,?/') 
Therefore, ([SD]) becomes 

Aij{Dij,^P)dijiP = Ti^P,Vij), (81) 

where 

J^{4^,V^J;) = F{^p)F'{lP)H{\V^P\^^|;) ( . (-^' + ^^'(Co - 1))|V^P ^2^ _ 
^ ^ ^' ' ^\^/72(|vV;|2,V;)c2-A(|VV|2,?/^) y 

Note that has quadratic growth in IVV'I, so it is not easy to get a prior estimate and 
the existence for (|8ip directly. The strategy here is that, instead of (|79p . we first solve the 
problem 

Aij{D^,7P)dij7p = T{iP,ViP) in Ql, 

X2-fl{xi) 
/2(a;i)-/i(xi) 



= iffi^"^ on af^i. 



(82) 



where 
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Thanks to (j7ip . one has 



(83) 



It follows from Theorem 12.5 and Remark in P308 in [T7] that there exists a solution tpi 
to ([82]) . Furthermore, writing -0^ = min{'0L,O} and t/ij = maxj-i/^^, 0}, by the proof of 
Theorem 3.7 in [T7\ . 



inf 'ipr — C sup 



< tpL ^ sup Tp^ + C sup 



where C = e'^ — 1 with d = sup{/2(xi) — /i(xi)}. Thus, 



_^^i-27 < < ^ + C<5^~2'^ for k sufficiently large. 



(84) 



Moreover, one can get some nice estimates for -0^. This follows from the techniques 
developed in Chapter 12 in [T7]. Using the specific form of estimate (12.14) in P299 in [T7] 
and Remark (4) on global estimate for quasiconformal mappings in P300 in ^17j, then one 
can improve the estimate in Line 7 in P304 in [17J to the following more precise form 



Mi,a<C(7,r?) l + |Z)u|o + 



(85) 



actually, C(7, fi) depends only on the diamfi and norm of 50. Here we use notations 
and symbols in (j85p as those in Chapter 12 in jl7j . 

Note that although the estimate (j85p is derived with zero boundary conditions, it holds 
in the case that the boundary value is constant in each connect component of boundary. 
Indeed, first, it holds for the case that the boundary value is a constant. Then one can 
generalize the estimate to the case that boundary value is constant in each connected 
component of the boundary, since all estimates are obtained through localization. 

Applying the estimate ([85]) to the problem ([82]) shows that, there exists = //(A/A) > 0, 
such that for any ^ and for ijjk with k > 4L, one has 



[V'fc]i,^;Bi(a.o)nQi < C{A/\ \fi\2) 1 + |^V'fclo;Bi(xO)nni, + 



(86) 
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This, together with interpolation inequahty and (|8ll) . yields 
<VC{A/X, \fi\2) 1 1 + \DiPk\o;Br{xO)nn^ + 





^ +Cn^m + C 










) 







where C, appeared in last term, is the same as that in (j84|) . Taking 770 sufficiently small so 
that 7/C(A/A, \ fi\2) < 1/2 if ?? < rjQ, then one has 



<r?C(A/A,|/,|2) + 
Thus, the Holder estimate (I86D becomes 



+ C;, m + C 



(87) 



'Jk\\i,ii,Bi{xO)nnL ^llV'fc||i;Bi(xf')nCi + [V'fc]i,fj;Bi(xO)nQL 

<(1 + C(A/A, |/.|2))||V'fc||i;Bi(.o)nn^ + C(A/A, l/.h) ( 1 + 



<C(A/A,|/i|2) (^r?oC(A/A,|/i|2) [ 1 + 

+ C(A/A,|/,|2) (^1 + 
<C(A/A,|/,|2)) ^l + m + 
Note that, for any x, y £ Ql, 

\vMx)-'^My)\ 



+ Cr,o\m + C 



< < 



\x - y[ 



ll^fe||i,M;Bi{x)nCi,, if y e ^1(2;) n VIl, 
2||'(Afc||i;ni, iiy i Bl{x)r^nL. 



This, together with ()87p and (j88p . yields the following Holder estimate 

m]i,f,;nL = sup ■ < C(A/A, \fi\2) 1 + m + 



Furthermore, it follows from ([88]), the Schauder estimate (Theorem 6.2 and Lemma 6.5 in 
|17j). and the bootstrap argument that 



'jk\\2,a;Bi,2{xO)nnL ^ ^ A/ X, \fi\c2.c , m, 
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Similar to the argument for (j89p . one has 

\\ipk\\2,a;nL < C A/A, \fi\c2,c,m, 



(90) 



Hence, using Arzela-Ascoh lemma and a diagonal procedure, we see that there exists a 
sequence ipk^ such that 

ipki ip in. C'^'^{K) for any compact set K <Z Cl and (3 < a. 

Furthermore, if) satisfies the problem 



3^2-/l(j^'l) 
/2(2:i)-/l{a;i) 



m on (90, 



+ C,, m + C 



and the estimate 

||V'l|i,n<r?C(7,|/.|2) (^1 + 

where 77 G (0,r/o). Thanks to estimate ([83]) . one has 

< r/C(7, |/i|2)(l + + C^(m + CJ^-^^), 

where C depends only on fh, A and A. 

Obviously, there exist tji G (0, r^o) and (5i € (0, 5o] such that 

ViCii, \fi\2)il + C5l~^^) < V(S2(^-eo)-2eo)/2, 
C^,(257/' + CSl-^^) < V(S2(fi-eo)-2£o)/2. 

Therefore, for any (5 G (0, and m S (5''', 2(5^''^), the solution ip satisfies 

<^\B-eo)-2eo. 



Now ([78]) follows from ([91]) and (l92]l . 

Furthermore, (|89|l and ([90]) yield the following higher order estimates 

\i,,;n<C{A/X,\M2) + m + 
\2,a;n < (^A/A, |/»|c2,«,m. 
This finishes the proof of the Proposition. 



(91) 



(92) 



(93) 
(94) 



□ 
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4 Far Fields Behavior, Existence and Uniqueness of Bound- 
ary Value Problem for the Stream Function 

In this section, we will study far fields behavior of the solution to (|77p. It will be shown 
that the flows induced by the solutions to ((771) satisfy asymptotic behavior (fT7I) - p^ . This 
also yields that solutions to (|77|) satisfy (|61|) . Combining (|61|) and (|92|) . we can remove 
both extension and truncation appeared in (j77p . Therefore, these solutions solve problem 
(j67p . Furthermore, the asymptotic behavior is crucial for our formulation for the problem 
since the stream function formulation is consistent with the original formulation of the 
problem for the Euler system in the infinitely long nozzle, as long as the flow induced by 
a solution to ([67|) satisfies (fT7|) - (fT8]) and ([28]) . Finally, the uniqueness of the solutions will 
be a consequence of the asymptotic behavior. To study the solution in its far fields, we will 
use a blow up argument and an energy estimate. 

For xi < n, define ip^'^\xi,X2) = Tp{xi-n,X2)x{fi{xi-n)<x2<f2{xi-n)}- For any compact 
set K (s (-00, oo) X (0, 1), it follows from ([MD that 

IIV'^"^ ||c2>«(i<') < C for n sufficiently large. 

Therefore, by Arzela-Ascoli lemma and a diagonal procedure, there exists a subsequence, 
ipi^k) ^ such that 

^ ^ in C^'^iK) (95) 

for any K d (— oo, oo) x (0, 1), for any /3 G (0, a). Furthermore, it follows from dH)-© and 
()94p that ^ = on X2 = and ijj = m on X2 = 1- So, satisfies 

^ ' {^^) = mmHim^ ^) m d, ^^^^ 

■0 = on X2 = 0, ■0 = m on X2 = 1, 
where D = (— oo, oo) x (0, 1). Moreover, by ([78]) . one has 

1-01 < C{eo, 6) and | V^l^ < ^^(5 - eo) - 2eo. (97) 

Thus, by the similar argument in Section [3l on any compact set E C (— oo, oo) x [0, 1], 

\mci..iE)<Cie,5). 
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Moreover, it follows from the Schauder estimate for second order uniformly elliptic equations 
that 

\mc2,c.^E)<C{e,5). (98) 
Therefore, G C'^'^{D). In fact, we have the following stronger results 

Lemma 5 There exists 82 G (0, 5q\ such that if 
(i)- ll^'llco.i([o,i]) = (5 < ^2, 

(a), m G {5'^,fh), where fh is defined in (56\) in Section\^ 

(Hi), there exists e < such that satisfies 

IV'I < C{e, 6) and | V^p - ^'^{Bii;)) < -e, (99) 

and solves the problem Ii96\) . where B is defined in Section\^ 

then tp is independent of xi, moreover, 

ip{xi,X2) = ■ip{x2) = poUo{s)ds, (100) 
Jo 

where po and uq are uniquely determined by B and m as in Section [B 

Proof: The proof is divided into two steps. First, it will be shown that %p is independent 
of xi. Then we will prove that ^ip is of explicit form (jlOOp . 

Step 1. Set w = ipxi- Differentiating the equation in ([96]) with respect to xi yields 



where Atj, 6 and ■(? are defined as 



(101) 



Ai,{q,z) = H{\q\^z)6ij-2Hi{\q\^,z)qiq„ (102) 
e(s,z) = {F"{z)F{z) + {F'{z)f)H{s,z)+F'{z)F{z)H2{s,z), (103) 
T?(s,z) = 2F{z)F'{z)Hi{s,z), (104) 



28 



for any g G s > 0, and z G M, here F" e L°°{M.^) since ||^'||co.i(R') - C6^~^^. It follows 
from (p9|) that there exists a constant A depending only on e such that 

\Ai,{Di;,i;)\<A{e). 

Although it is unknown whether -0 G C^(D), the equation (jlOip holds in weak sense. 
Moreover, ij; satisfies the boundary conditions 

= on X2 = 0, 1. 

Let ry be a Cq" function satisfying 

r] = 1 ioT \s\ < I, rj = for \s\ > I + 1, and |r/'(s)| < 2. (105) 
Now multiplying rf'{xi)w on both sides of pUl|) and integrating it over D yield 



0(1 V'(/'p, V')??^if^c?a;ida;2 + ^{\V4>\'^ ,4')di'il^diwr]'^wdxidx2- 

D 

Substituting the explicit forms of Aij, i^idV-i/^p, -i/i) and i^2(|V'(/5p, -(/J) into the above equal- 
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ity and noting that if; satisfies (|99|) . one may get 



DHi\V^\^,^) 



dxidx2 



JJd H^{\V^\Mj) ' ^ 1/12 
2— —OjWOiri'nwaxiax2 



F'{i))F{4))H2{\Vi}\'^,i})ri^w'^dxidx2 
2F'(Vi)F(Vi)#l(|VV^|^V')^^^fV^ • Vwdxidx2 

dxidx2 



2\ 



D 

D ^(|V'0|2,V;)(^2(|V^|2,^)c2 _ IVV^I 

2__«i( — ^'^j^ djwdirirjwdxidx2 
D H'^{\V%1)\'^,%1)) 

yyD//2(|VV;|2,V)c2-|V^|2 ' 



2 



D i72(|W|2,V)c2 - 



which can be written as 

iff\Vw\ 



-dxidx2 = ^^-^j- (106) 



First, it is easy to see that Ii + + < 0. Second, due to (ITTI) . one has 



l/sl < C5^-^^ I / w^dxidx2. (107) 



JO 
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Finally, since H < g{B), thus if 62 is sufficiently small, one gets from above that 
ri n 



\Vw\ dx2 



/ dxi / 

J -I Jo 

< C{B,e)( I dxi+ I dx^ [ \Vw\'^ + \Vww\ + w'^dx2 

\J-i-i Ji J Jo 

+C{B)6^-^^ [ [ w^dxidx2 
J -I Jo 

< C{B,e)( f dxi + dxi^ [ \Vw\'^ + w'^dx2 + C{B)6^-^^ [ [ w^dxidx2. 

\J-l-l J I J Jo J -I Jo 



Notice that if = on 2:2 = 0. It follows from Poincare inequality that 



1 ,-\ 

v?dx2 < / \Vw\'^dx2. (108) 



'0 Jo 

Therefore, there exists a constant C independent of / such that 



f I \Vw\'^dxidx2 <c( I dxi+ f dx^ I \Vw\'^dx2 (109) 
J -I Jo \J-l-l J I J Jo 



for large /. It follows from (j98p that 

-I rl+l \ rl 

dxi + dxi] \Vw\'^dx2 < C 
Ji J Jo 

for some uniform constant C independent of /. Therefore, 

dxi / \'Vw\'^dx2 < C 
-I Jo 

for some constant C. Taking / — > 00 yields 

dxi / \'Vw\'^dx2 < C. 
-00 Jo 

Hence 



/dxi + / dxi j / \\7w\'^dx2 ^ as Z — > cxd. (HO) 
Ji ) Jo 



Taking the limit Z — > 00 in (jlOOp . one has 

r"oo c\ 



00 Jo 



\\Jw'^dx\dx2 = 0. 



So 

w = 0. 
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(Ill) 



Therefore, ip = ^{x2). Thus ip solves the following boundary value problem 

ipiO) = 0, = m. 

Step 2. Uniqueness of the solution to the boundary value problem (lllip . 
Suppose that there are two solutions 'tpi and 'tp2 to (|llip . Let = — ■02- Then 
satisfies 

{a(j)' + 60)' = c4>' + d4>, 
m = 'A(l) = 0, 



112) 



where 







-ds, b 



H^m,i:) 



1 - 

^{\il>'\'^,il>')^'ds, d= / Q{\il,'\'^,^')ds, 
Jo 



with ip = sipi + (1 — s)^2; where and are defined in (jl03p and (jl04p respectively. 
Multiplying on both sides of the equation in (jll2p . and integrating it over [0, 1], we have 



Jfj' dx2<- /V(F^W)^ + FWF^^W)ff(|V^f,^)<^^dx2. 
H{\'ip'\'^,ip) Jo ^ ' 

Note that ||i*"||cO'i(Ri) — C8^~^'^ , thanks to the smallness of b and the Poincare inequality, 
one has 



< 0. 



Therefore, </> = 0. So the solution to (jllip is unique. On the other hand, by the definition 
of H and -F, one knows that the boundary value problem (jllip has a solution 



iP = ^(X2) = / PQUQ{s)ds. 

Jo 



This finishes the proof of the Lemma. 



□ 



It follows from Lemma [5] and (j95p that the fiow induced by the stream function satisfies 
([HD and (HI]). 

The asymptotic behavior in the downstream can be obtained by a similar argument. 
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An important direct consequence of this far fields behavior is a better maximum estimate 
for the stream function. 

Proposition 6 // ||5'||co>i([o,i]) = S < mm{6i,62}, B'{0) < 0, B'{1) > and m e {6^,m), 
then the solution established in Proposition^ satisfies (E2P- 

Proof: It follows from Proposition [5] that 

ip{xi,X2) — > / poUo{s)ds uniformly as xi — oo, 
Jo 

'0(xi,X2) I piui{s)ds uniformly as xi +oo. 

J a 

Therefore, for any e > 0, there exists L > such that 

— e < tjj{xi,X2) < m + e if > L. (113) 
Note that F'{i{j) > in the domain {ip > m}, thus 

Aij{D'ilj,'ip)dij'ip > 0, in the domain {ip > m} D {\xi\ < L}, 
where Aij is defined in (j76p . By maximum principle, one has 

— e < ip{xi,X2) < m + e in {ip > m} n {\xi\ < L}. 
Since F'{%1)) < in the domain {tp < 0}, thus, similarly, one can show that 

— e < ij{xi,X2) < m + e in {^/i < 0} n {|xi| < L}. 
Combining these estimates with pi3p . it yields 

— e < ip{xi,X2) < m + e in ri. 

Since e is arbitrary, one has 

< ^jJ{xl,X2) < m in 0. 
This finishes the proof of the Proposition. □ 
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It follows from estimates ([6T]) and ([78]) that the solutions established in Proposition d] 
are solutions to ([67|) when the assumptions of Proposition [6] are satisfied. 

In fact, one can also use energy estimates to show that uniformly subsonic solution to 
([67j) is unique. 



Proposition 7 Let the boundary of 0, satisfy Then there exists 5^ £ {0,5q] such 

that if 

(i) - ll^'llco.i{[o,i]) = ^ <h, 

(ii) . m G {5"' ,7n), 

then there exists at most one solution %[) to ( [g7[ ) satisfying 

< V < m, IVV'P - t^iBii))) < -e for some e > 0, (114) 
where H and F are defined by B and m as in Section\^ and B{tp) = /i(yOo) + — 
Proof: Let tl^i and ip2 be two solutions to (j67p . Set ip = ipi — ■iIj2- Then -0 satisfies 



di{aijdj'il)) + di{bi'4)) = Cidiip + d^, 
= on U ^2, 



(115) 



where 



H^{\V^P\,^P) ' ' Jo H^{\V^P\,^P) 

1 rl 

^{\Vi|Jf,i;)^^i;ds, d= / Q{\Vil,\\i,)ds, 

JO 

here ■0 = s-^i + (1 - s)'02, -4jj, 6 and ^ are defined similar to piI2|) . pi]5|) and piH]) . 
respectively except we replace F and H hy F and ff. 

Multiplying r]'^^p^ on both sides of equation in (jll5p . where r] is defined in ()105p and 
'(/'^(x) = max{'0(3;), 0}, then similar to the proof of LemmaO one has 

[[ \Vil)\^dxidx2<C{B,e) f f \Vi)\^dxidx2. 
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It follows from Lemma [5] that ipi and ^p2 have the same far fields behavior, thus and 
IVV'I ^ as |xi| ^ oo. Thus 



inn{ip>o} 

so ^ 0. Similarly, one can show that ^ > 0. Therefore, ^ = 0. This finishes the proof of 
the Proposition. □ 



5 Refined Estimates for the Boundary Value Problem for 
Stream Functions 



In this section, we will derive some refined estimates for solutions to the problem (I67p . 
Combining these refined estimates with the estimates obtained in section [3] and section [U 
one will get a solution to the Euler equations ([I])- ([3]), with the boundary condition ([8]) and 
the constrains ([9]) and ()13p . More precisely, it will be shown that ^lJx2 is always positive, 
therefore, u = ipx2/p = V'x2/-f^(l ^V'Pi "0) satisfies p^ . This positivity of the horizontal 
velocity and the asymptotic behavior yield that (p, n, v) induced by ip satisfies the original 
Euler equations, the boundary conditions and the constrains on mass flux and Bernoulli's 
constant. 

Lemma 8 Let the boundary of 0, satisfy Then there exists 84 G (0, Jq] such that if 

(i)- ll^'llco.i{[o,i]) = ^ 
(a), m £ (5'^ ,fh), 

(Hi), ip satisfies 1^114^ and solves {61) with H and F defined by B and m as in Section\^ 



then V' satisfies 

< ip < m in Q, (116) 

and 

1^X2 >0 inCl. (117) 
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Proof: It follows from (jll4p that -0 achieves its minimum on Si and maximum on 52, hence 

V'xa > on on. (118) 

On the other hand, U = ipx2 satisfies 



(119) 



= e{\Vtp\'^,i;)u + mvi^\'^,i^)di^diU 

in the weak sense, where Aij, Q and •& are defined similar to (|102p . (jl03|) and (|104p respec- 
tively except we replace F and H hy F and We first claim that 

[/ > in Q. (120) 

Indeed, it follows from Lemma[5]that U{xi,X2) > when > I for some / sufficiently 
large. Multiplying (11190 by = min{C/, 0}, and using (IllSp . one may get that 

e(|VV'P,V')t/^(ixidx2 - // i}{\Vi}\'^,i;)UVij ■VUdxidx2 
{u<o} JJd 

< -II {F"{^)F{il,) + {F'{^)f)H{\Vil>\\^lj)U^dxidx2 

J J{U<0} 

< C5^-^^ 1 1 U^dxidx2. 

J J{U<Q} 

Define Kx^ = {2;2|/i(a;i) < < f2{xi),U{xi,X2) < 0}, then K^^ is an open set for each 
xi. Let = UiQ.yll^, where I*^ are connected components of K^^^- For each X2 G 

U{xi,X2) = / U{xi,s)ds. 

J min /i. 
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Therefore, 



{U<0} 



U'^{xi,X2)dxidx2 



/ dXiT,i^y / U^{xi,X2)dX2 

I dxiT,i(zy I i I dx2U{xi,s)ds] dx2 



1 \ """'^i 
I r />max/!. 



< / dxiSjg^' / / {dx2U{xi,s)) (is(max/*^ — mini* Jci3;2 

fl /-max It 



dxiSjg.y(max/*^ — mini* / {dx2U{xi, s))'^ds 

-I J rain 11^ 

dxiY,i(zy / {dx2U{xi,s))'^ds 

-I J rain 11^ 

< max |/2(xi) - [ [ \VU\^dxidx2 



Hence, 



{U<0} 



dxidx2 < 06^-'"^^ 1 1 \VU\^dxidx2, 



{c/<o} ^^(|VV'P,V') JJ{u<o} 

which imphes 

/ / \VU\^dxidx2 < 0, 
J J{U<0} 



so must hold. 

Now, we use an argument similar to the proof of Lemma 1 in §9.5.2 in ^llj to show that 

= [/ > in O (121) 

holds for any weak solutions U to (jll9p . 

Indeed, let U = e~^^^U. Then f7 is a nonnegative weak solution to 

di (^^e'^^^d.U^ + - - ^(1 VV'|^ e-^a,C/ + Ge^^C/ = 0, 

where Aij and i) are defined in (|102p and (jl04p . and 

^22 2 , fr. f A2\ H2d2tp A „ (H2di'4)\ , ,2 
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with G defined in (jl03p . Choose cr > sufficiently large. Then G > 0. Thus 

^<^^2f,T~T\ , Mi2^ ^25iV ..|V7„/,|2 



It follows from Theorem 8.19 in [17] that (fml) holds. 
Now, (fTT6|) follows directly from (fTTi]) and (fT2T]) . 

Since ip = m on S2, if F'{m) > 0, then for any {x\, f2{x\)) G £'25 there exists a small 
disk c r2 satisfying J\f r\Q. = {x\, f2{x\)) such that F'{tp) > in 7\A, therefore, 

Aij{Di},'4))diji) > in AA. 

Moreover, by (|116p . < m in J\f . Thus, by the Hopf Lemma, one has ipx2{x\-, f2{x\)) > 0. 
In the case F'[m) = 0, ip satisfies 

Aij{Dip,il;)dij{ip - m) + R{ijj - m) = 0, 

where 

F{'p)H^c^ F'{ip)-F'{m) 
~ ~i72(|V?/^|2,V^)c2 - IWP i^-m ■ 

It follows from the Hopf lemma that 

dx2'4' > in ^2. 

Similarly, one can show that 'ipx2{xiy fii^i)) > for any xi € M. 

This finishes the proof of the Lemma. □ 



Choose 5o = m[n{5i,62,6s,64}, then 5o > 0. If ||-B'||(^o,i([o,i]) = ^ < <^0) for any m E 
{6^ , 26q^'^), there exists a solution to the problem (j67p . It follows from Lemma [5] and Lemma 
[8] that the flow fleld induced by tp satisfies (j28|) and ([29|) . hence Proposition [2] guarantees the 
existence of Euler flows. Furthermore, Proposition [2] and Proposition [7| imply uniqueness of 
Euler flows with asymptotic condition (jlSp . mass flux condition (jlSp . and asymptotic 
behavior (fT7l)-(fT8]). 
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6 Existence of Critical Mass Flux 



So far, we have shown that, for the given Bernoulh's function in the upstream satisfying 
(jl4p . there exist Euler flows as long as m G {S"^ ,25^^"^). In this section, we will increase m 



as large as possible. 

Proposition 9 Let $7 satisfy and B satisfy and [21\) . Then there exists rh <fh 

such that if m ^ {5^ ,m), there exists a unique ip which satisfies 

{) < < m in n, and M{m) = sup {\Vi;f - ^^(^(V'))) < 0, (122) 

n 

and solves (67^, where B{ip) = h{po)+F'^{Tp)/2. Furthermore, either M{m) as m ^ rh, 



or there does not exist a > such that {67} has solutions for all m G {rh,m + a) and 



sup M{m) < 0. (123) 

mg(rTi,m,+(T) 

Proof: The basic idea of the proof for Proposition is quite similar to that in [21 [31] . 

For the given Bernoulli's function B in the upstream satisfying (I14p and any m G (5^ ,fh), 
one can define po and uo{x2), and therefore F{'ip) as in Section [2l Note that po s^nd F 
depend on m by definition, thus in this section, we will denote them by po{m) and F{^p■, m) 
respectively. 

When B satisfies (I21D. one has 



F'{m) = F'(0) = 0. 
Thus F', the extension of F' is Section 3, has the following simple form 



F'{s) 



F'{s), ifO<s<m, 
0, if s < or s > m. 



(124) 



Set F{s) = F'{s)ds. Then it is easy to check that 



Bo<B< h{po) + ^^Y^ < B and ||F'||co,i(mi) < C^^'^T'. (125) 
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Let {en}^=i be a strictly decreasing sequence of positive numbers such that ei < eo/4 
and £n i 0. One can truncate H associated with e„ as follows 

H^^\\V^\^t,m) = J{An{\ViP\\tl;;m),Bn{ip;m)). (126) 

To give a clear explanation of this definition, we first introduce a sequence of smooth 
increasing functions Cn such that 



Cn(s) = < 



s, if s < —2en, 

-£n, if S> -£n- 



Now one can define 
where 

Bn{W]m) = h{po{m)) H . 

It is easy to see that there exist two positive constants A(n) and A(n) such that 

A(n)|e|2<i;;)(g,z;m)e.e, <A(n)|e|2 

for any z G R-*^, g G and ^ G R^, where 

A^\q,z;m)=H^^\\q\^,z;m)dij-2H[''\\q\'',z;m)qiq,. 

Thus it follows from the argument in Section [3] that for any m G ((5'^,m), there exists a 
solution ■0("')(x;m) to the problem 



A^^\Dil),'il);m)dijilj = Tn{DTp,ip;m) in Q, 



(127) 



where 

where we ignore some obvious independent variables in the definition of J^n- Moreover, if 

|VV^(")|2 -^(^(");m) < -2en, (128) 
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then = 1. Similar to Section [3l one has 

< ^'^"■^(x;m) < m. 

Since F satisfies (jl25p independent of En, one can estimate I5 in (|106p as that in (jl07p . 
Furthermore, it fohows from the same arguments in Lemma [5] that the solution to (jl27p 
satisfying (1128P has far fields behavior as (jlOOp . In addition, by Proposition [7l such a 
solution is unique among the class of solutions satisfying pOUp . 

Note that in general, we do not know uniqueness of solutions to problem ()127p . Set 

5„(m) = {V'^"^(x;m)|V'^"n2;;m') solves the problem (fT27l) }. (129) 

Define 

Mn{m)= inf sup(|VV'(")(x;m)|2 -S2(^(V'(");m))), (130) 

and 

Tn = {s\6^ < s < rh,Mn{m) < -4e„ if m G (5^,s)}. 

It follows from Proposition HJ Lemma [5] and Proposition [6] that [(5q,2(5J^^] C Tn, therefore, 
Tn is not an empty set. Define m„ = supT„. 
The sequence {m„} has some nice properties. 

First, M„(m) is left continuous for m G ((5'>',m„]. Indeed, let {ml^^} C {6'^,mn) and 
mi^^ t m. Since M„(m^'^^) < — 4e„, one has 

||^(")(x;mW)||^...(n)<a 

Therefore, there exists a subsequence ■ip^'^\x;m^n'^) such that 'ip^'^\x;ml^'^) — > tp, moreover, 
solves (fT271) . Thus M„(m) < limM„(mi'''^). So Af„(m) < -4e„. Note that all these 
solutions satisfy the far fields behavior as (jlOOp . by uniqueness of solutions in this class, 
M„(m) =limAf„(m^f)). 

Second, m„ < fh. Suppose not, by the definition of m„, fh S T„. It follows from the 
left continuity of Mn, Mn{fn) < — 4e„. Thus by means of the proof of LemmaO '4>^"'\x;fh) 
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has far field behavior as in (jfOOp . However, it follows from the definition of m that 
sup (|VV'(")(x;m)|2 - s2(i3„(V'("^(x; m)))) 
> sup m.aii{{\pQ{fh)uo{x2;rh)\^ - i:^{B{x2))), {\pi{m)ui{y{x2);m)\^ - T.^{B{x2)))} 

2^2G[0,l]) 

= 0, 

where y = y{s) is the function defined in (j49p . Thus Mn{rn) > 0. This is a contradiction. 
Therefore m„ < m. 

Finally, {mn} is an increasing sequence. This follows from the definition of {m„} directly. 
Define rh = lim„^oo"i-n- Based on previous properties of rh is well-defined and 

rh <rh. 

Note that for any m G (5'>',m), there exists > m, therefore Mn{m) < — 4e„. Thus 
ip = ■(/;("■) (x; m) solves (j67|) and 

supdVVp - = M„(m) < -4e,. 

If sup^g(57 M{m) < 0, then there exists n such that sup^g(-57 M{m) < — 4e„,. As 
the same as the proof for the left continuity of Mn{m) on (5''', m„], Mn{rri) < —Ae. Suppose 
that there exists a > such that (j67|) always has a solution for m G (m, m + o"), and 

sup M(m) = sup supdVV'P - < 0. (131) 

mG(m,m+CT) mG(m,rn+CT) f2 

Then there exists > such that 

sup M{m) = sup supdVV-l^ - s2(i3(V'))) < -4e„+fc. 

m£(m,m+(T) mS{m,m+(7) H 

This yields that m„+fc > m + a. So there is a contradiction. The contradiction implies 
that either M{m) — > 0, or there does not exist cr > such that (I67p has solution for all 
m G (m, m + fj) and (I123P holds. 

This finishes the proof of the Proposition. □ 

It follows from Lemma [5l Lemma [8] and Proposition [9] that if B satisfies ()14p and 
m G ((5''',m), then there exists an Euler flow through the nozzle. Collecting all results 
obtained together, we complete the proof of Theorem [TJ 
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